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Abstract 

Certain regions of the solar atmosphere are at sufficiently 
low temperatures to be only partially ionised. The lower 
chromosphere contains neutral atoms, the existence of 
which greatly increases the efficiency of the damping 
of waves due to collisional friction momentum transfer. 
More specifically the Cowling conductivity can be up to 
12 orders of magnitude smaller than the Spitzer value, 
so that the main damping mechanism in this region is 
due to the collisions between neutrals and positive ions 
(FKh odachenkoet al. 200 41). Using values for the gas 
density and temperature as functions of height taken from 
the VAL C model of the quiet Sun (IVernaz za et al 19811 ). 
an estimate is made for the dependance of the Cowling 
conductivity on height and strength of magnetic field. Us- 
ing both analytic and numerical approaches the passage 
of Alfven waves over a wide spectrum through this par- 
tially ionised region is investigated. Estimates of the ef- 
ficiency of this region in the damping of Alfven waves 
are made and compared for both approaches. We find 
that Alfven waves with frequencies above 0.6 Hz are 
completely damped and frequencies below 0.01 Hz un- 
affected. 



1 Introduction 

MHD waves are vital phenomena in the study of the 
solar atmosphere. They have been suggested as a 
source of heating of both the chromosphere and the 
corona ( IGoodman 20001 . |Piddington 195 6^), based on 
the damping of wave energy due to dissipation effects. 

Three main energy dissipation mechanisms in the chro- 
mosphere have been studied in previous work, namely 
collisional effects, viscous effects and thermal conductiv- 
ity effects. Viscous damping of MHD waves is caused by 
the momentum transfer from the thermal motion of par- 
ticles. Collisional friction is caused by the relative ve- 
locities of the species in the plasma. While thermal con- 
ductivity is also related to kinetic effects and momentum 
transfer Khodachenko et al. (2004) and Khodachenko et 
al. (2005) showed that in the partially ionised plasma of 
the lower solar chromosphere, collisional effects are far 
more important than the other two. Also the presence of 
neutrals means that of all the collisions to be considered 
it is the collisions between neutrals and ions that are the 
most important component of this mechanism. 

Alfven waves are a specific type of MHD wave. 
Previous work on the generation and propagation 
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of Alfven waves has included the simulation of 
high frequency waves generated by photospheric mo- 
tions and the damping of these waves in the par- 
tially ionised chromosphere above. This damp- 
ing was suggested as an explanation of the for- 
mation of spicules (IDePontieu and Haerendel 19981. 
PePontieu 19991, | James et al. 2003 1). Spicules are 
high, thin jet like structures consisting of chromospheric 
plasma. They can reach heights of between 5000 and 
9500 km whereupon they either fall back or fade in the 
hot corona. This work involved the numerical simulation 
of Alfven waves of frequencies around 0.5 Hz damped by 
ion-neutral collisions in the WKB approximation. This 
approximation assumes that over a single wavelength 
there is no spatial variation of variables. More specifi- 
cally, the wavelength itself does not change much on its 
own scale and the damping scale is much larger than the 
wavelength. De Pontieu et al. (2001) analytically calcu- 
lated the damping time for Alfven waves of varying fre- 
quencies in model chromospheres of various solar struc- 
tures, given estimates for the plasma parameters in these 
regions. They estimated that frequencies above about 0.1 
Hz are unable to penetrate through to the corona from the 
photosphere. 

Other types of MHD waves in the solar atmosphere 
have also been studied. Linear analysis of slow MHD 
waves in the frequency range below 3 .5 x 1 0"^ Hz has been 
used to show the damping as a possible heating mecha- 
nism for the chromospheric network. The dissipation of 
the currents asscoiated with these waves, which is sug- 
gested to occur in thin magnetic flux tubes with strengths 
of 700G to 1500G, heats the tubes and hence the chromo- 
spheric network (| Goodman 2000 1). 

The passage of waves through both the chromosphere 
and corona has been studied using measurements of 
magnetic bright point (MBP) positions in the photo- 
sphere. A frequency power spectrum for horizontal mo- 
tions ( |Cranmer and Ballegooijen in pre ss)) at the photo- 
spheric level were derived in the frequency range 10"^ 
to 0. 1 Hz. Using these as a lower boundary condition, a 
WKB approximation was applied to derive power spectra 
at different heights in the atmosphere to show the effective 
damping of waves in different regions. It was shown that 
waves in this range (10^^ to 0.1 Hz) may be evident up 
to a few solar radii which suggests that horizontal pertur- 
bations in this frequency range may be unaffected by the 



partially ionised chromosphere. 

Higher frequency waves have also been detected in 
the solar atmosphere. Waves of frequencies up to 
0.1 Hz have been detected in the upper chromosphere 
( IIDefores t 20041) and transition region. These waves 
were visible in the TRACE 1600 A passband. The ex- 
act source of these waves is uncertain and they are not 
energetically significant to the chromosphere. High fre- 
quency waves (0. 1 Hz) in the lower chromosphere have 
proved harder to find ( HFossum and Carlsson, 2005) ). It 
has been suggested that high frequency Alfven waves in 
the range 1 to 800 Hz can be created by small scale recon- 
nection in the chromospheric network, and their dissipa- 
tion was considered as a possible heating mechanism for 
the corona above ( yviarsch and Tu, 1997j ). 

Here we present both analytical estimates and full 
MHD simulations of the efficiency of damping of Alfven 
waves in the partially ionised chromosphere. We give es- 
timates of the effective filter function for the damping of 
Alfven waves. This damping is assumed to be due to the 
collisions of neutrals and positive ions in a model hydro- 
gen plasma. The Alfvenic disturbances are assumed to be 
created at the photospheric level by horizontal motions on 
the surface. 

Given values for the plasma density and temperature in 
the lower chromosphere taken from the VAL C model of 
the quiet Sun (TV ernazza et al 19811 1 the Cowling con- 
ductivity is estimated. Using the hnear damping approx- 
imation (IBragins kii 1965| ) the damping decrement for 
upwardly travelling waves is calculated and the efficiency 
of this region in damping the waves is estimated. 

Then full MHD numerical simulations of linear Alfven 
waves in a model atmosphere which contains a partially 
ionised chromosphere are performed. We investigate the 
efficiency of collisional damping for a range of frequency 
of wave. The numerical results are compared with the 
theoretical estimates and are shown to be in agreement 
with the linear damping approximation for low frequency 
waves but differ for higher frequency waves. 
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2 Collisional friction damping in a 
partially ionised plasma 

To incorporate the damping of MHD waves due to ion- 
neutral collisions in a partially ionised plasma, the gen- 
eralized Ohm's law should be included in the governing 
equations. ( |Cowhng 1957| , fBraginsk ii 1965| ). 
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where v is the fluid velocity, B is the magnetic field,E is 
the electric field, j = (V A B)/jU() is the current density and 
jj. is the component perpendicular to the magnetic field. 
The conductivity is defined by 
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The number densities of species (ion, electron, neutral) 
are given by ni,ne,n„ respectively, and the masses by 
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is the relative density of neutrals, and v',v, v',v and v',(, are 
the effective collisional frequencies defined by 
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where k = e,i,l - i, n. Following the example of 
fSpitz er 1962| , the collisional frequencies of ions and 
electrons with neutrals are estimated by 
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and Kb is Boltzmann's constant. The collisional fre 
quency of electrons and ions is given by 
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( ISpitzer 1962| ). Z,-,, = 5.10-iW,S,„ = lO-^'^m^ are 
the ion-neutral and electron-neutral cross-sections respec- 
tively. Here A is the Coulomb logarithm. 

In the generalised Ohm's law (equations 1 and 2) we 
have neglected the pressure term, as the chromospheric 
plasma is relatively cold, and the Hall term. The Hall 
term can be dropped from the generalised Ohm's law if 
the plasma is magnetised, i.e. if the ions and electrons are 
tightly bound to the magnetic field. In terms of the ion- 
gyrofrequency and the collision time this condition can be 
written as 

W,T»1. (11) 

Using the definition of the ion-neutral colUsion frequency 
this is equivilant to 
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l6kBT n„I,„ 

which depends on the strength of the magnetic field as 
well as the temperature and number density. Using val- 
ues for temperature and number density taken from the 
VALC model of the Sun ( IVemazza et al 1981 1), and var- 
ious models of magnetic field strengths in the chromo- 
sphere (see figure 1), this condition is satisfied every- 
where in the chromosphere. Khodachenko and Zaitsev 
(2001) showed that this condition does not hold in the up- 
per photosphere where the plasma is more weakly ionised. 
In section 3 we will show that the damping mechanism we 
are considering is not effective in the photosphere, and 
thus we only study Alfven waves in the chromosphere for 
which the neglect of the Hall term is valid. 

If the plasma is assumed to be entirely composed of 
hydrogen then the expressions can be simplified by taking 
min - m„/2, ^„ - pn/p resulting in 
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Using the definition for the Cowling conductivity as 
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and defining the Coulomb and Cowling resistivity as r] ■ 
l/cr and rj^ - 1/cr^, respectively it is trivial to show that 



Thus equation 2 can be be written as 



(15) 



E + (V A B) = m + (j]c- n)j± 

= m + VcL (16) 

and the frictional Joule heating term for this plasma 
( |Braginskii 1965J ) is then given by 
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with jx and jy being the components of j perpendicular 
and parallel to B. Thus in a partially ionised plasma the 
Coulomb resistivity acts to diffuse currents parallel to the 
magnetic field and the Cowling resistivity acts to diffuse 
currents perpendicular to the magnetic field. The damp- 
ing decrement (6) for this effect was derived by Bragin- 
skii (1965) for Alfven waves propagating with angular 
frequency ll) and wavenumber k and is given by 
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Khodachenko et al. (2004, 2005) showed that for typical 
magnetic field strengths and using the VAL C model for 
the plasma parameters, the ratio of Coulomb to Cowling 
resistivity can vary from 1 at the solar surface (z=0) down 
to 10 at heights of 2000 km. For the partially ionised 
chromosphere 
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Hence when considering frictional damping mechanisms, 
the ion-neutral colhsions are by far the most important. 

3 Analytic approach 

In order to evaluate the expression for the COwling con- 
ductivity an estimate for the neutral fraction ^„ is re- 
quired (see equation 15). Although the VAL C model 



gives us the ionisation degree for the Sun our model is a 
hydrogen plasma. Thus ^„ is calculated from the density 
and temperature values only of the VAL C model. Here 
we follow the method of De Pontieu (1999), assuming a 
hydrogen plasma where the number of electrons is equal 
to the number of protons (n, = n^,). 

The solar chromosphere, and spicules especially are 
not in LTE. Hence a simple one-level model for 
the hydrogen atom is inadequate for these conditions 
([Pottasch and Thomas 1959 1). A two-level model is used 
instead for the hydrogen atom, so that the ionisation equa- 
tion ([Brown 1973 1) is 
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where n is local number density, and x is the ionisation 
degree. Cj and Pj are the ionisation rates from level j 
due to thermal collisions and local radiation field respec- 
tively. Cj* and Pj* are the ionisation rates by non-thermal 
particles and external radiation field. In the first approxi- 
mation photoionisation from level 2 to level 1 is provided 
by the external field alone, so that P\ - P\* - 0. This 
external field is simply the photospheric radiation field at 
temperature Tr, suitably diluted by factor w. The ionisa- 
tion equation can be further simplified by noting that ther- 
mal collision ionisation is unimportant when compared to 
photoionisation (| Ambartsuyan 1958 1). Thus the ionisa- 
tion equation is balanced by photo-ionisation from level 
2 to level 1 and spontaneous recombination for the return 
route. 

dx 
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The steady state solution to this equation is given by 
( IThomas and Athay 196T| ) 
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where Tr is the temperature of the photospheric radiation 
field and w = 0.5 is its dilution factor. 

Using the solar atmospheric plasma variables given by 
the VAL C model the neutral fraction ^„ and the profile of 
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T/c with height can be calculated for a given magnetic field 
profile. 

As we intend to simulate the motion of waves in the 
chromospheric region of the solar atmosphere, we choose 
as our magnetic field model a spreading vertical flux tube 
representing an open magnetic structure with its footpoint 
in the photosphere. As Alfven waves propagate along a 
fieldline we can restrict ourselves to one spatial dimen- 
sion and need only define the variation of the magnetic 
field as a function of distance along the fieldline. Thus the 
horizontal expansion of the tube does not enter explicitly 
into our calculations. We choose a magnetic field profile 
to reflect the conservation of flux as the tube expands. A 
power law dependance on density is chosen 
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where pph is the photospheric density 2.7 x 10""^kg/m^ 
and Bo is 1200G. This profile is chosen so that at the 
base of the domain the magnetic field is approximately 
lOOOG and falls to a value of about lOG at a height of 
3000km. This is a simplistic model for the magnetic field 
but captures the decrease with height and enables compar- 
ison of analytic and numerical results. Martinez-Pillet et 
al (1997) showed that the field strength in flux tubes can 
drop to as little as 300G at heights of only 300km so we 
use this as the minimum value of magnetic field at this 
height for our models. The magnetic field profile for four 
diff'erent power law profiles are shown in figure 1 . A typ- 
ical calulated value of rfc as a function of height for the 
case when a - 0.3 is shown in figure 2. Khodachenko 
et al. (2005) calculated ^ for the solar atmosphere given 
constant magnetic fields. The ratio at a height of 1500 
km can be as low as 10"^ giving a value of the Coulomb 
resistivity of 8 x 10"^m^/s. 

Equation 18 gives the inverse damping timescale of 
Alfven waves with angular frequency lo and wavenumber 
k\\ along the magnetic field, propagating through a region 
of constant rjc- 

i = 2w5 = —ki (26) 



Figure 1: The variation of the model background mag- 
netic field with height. The four curves relate to the value 
of the exponent a - 0.3 (small dash line), 0.35 (big dash 
line), 0.4 (dot-dash line) and 0.6 (solid line). 
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Figure 2: Cowling resistivity as a function of height cal- 
culated from VAL C model of the quiet Sun and varying 
Bo with height according to equation 23 with a = 0.3. 



The linear MHD analysis used to obtain damping 
decrements (| |Braginskii 1965| ) assumes that the damping 
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decrement 6 is much less than 1. This impUes that 
S = « 1 



Damping Efficiency 



or that 
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For the model used here both ijc and the Alfven speed 
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are functions of height so the damping decrement also 
varies with height. Over a small timestep dt the amph- 
tude of the propagating wave changes by 
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By integrating over the partially ionised region where rjc is 
non-zero (figure 2) we get an estimate for the total change 
in amplitude due to the damping mechanism 
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where the integral can be performed numerically. The re- 
sults from this approximation are shown in figure 3 for the 
four diff'erent magnetic field models (different a). From 
these estimates it appears that in the absence of any other 
damping mechanisms and stratification efl'ects, high fre- 
quency Alfven waves are unable to pass through this par- 
tially ionised region without being completely damped. 
For the magnetic field model with a = 0.3 this can be as 
low as 0.1 Hz. 



4 Numerical method 

The propagation of Alfven waves in a partially ionised 
plasma is modelled in ID by numerically solving the non- 
ideal MHD equations, given here in dimensionless form. 



dp 



(32) 
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Figure 3: Analytical estimates for the damping efficiency 
of the partially ionised region of the solar chromosphere 
as a function of frequency. The four lines represent the 
four magnetic field profile given by a = 0.3 (solid line). 
0.35 (dash fine), 0.4 (dot-dash fine) and 0.6 (triple dot- 
dash line). 
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which also include the Coulomb and Cowling resistivities 
J] - Ijcr and rjc - I jcTc. B is the magnetic field, j is the 
current density, v is the velocity, P is the thermal pres- 
sure, and e = Pjiy - l)p is the specific energy density. 
The constants no and y have the standard meanings of 
magnetic permittivity of free space and the ratio of spe- 
cific heats respectively. Although for a partially ionised 
plasma y < 5/3 we found the value had little effect on the 
results so 5/3 was used in all simulations. 

It is worth noting that no viscous effects are included 
in these model equations as only coUisional damping is 
being investigated in this work. 

The simulations are carried out using a 2D MHD shock 
capturing code (Lare2d), applied in ID only. It uses a 
staggered grid with density, pressure and specific energy 
density defined in the cell centres, magnetic field at the 
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Figure 4: The background plasma temperature (solid 
Une), density(dot-dash) and magnetic field (dash), where 
a = 0.3. The values are all given in dimensionless units. 



cell faces and velocities at the vertices. For further de- 
tails see Arber et al. (2001). The code is adapted so that 
the resistive term in the induction equation updates the 
perpendicular and parallel components of current density 
separately, according to equation 16. Although the code 
is ID all vector variables have 3 components. 

The simulation domain extends 3000 km vertically, 
while the computational grid consists of 2000 cells. The 
number of cells vertically is restricted by the need to re- 
solve the smallest wavelengths in the parametric study. 
The vertical boundary conditions are line-tied which are 
perfectly reflecting. Thus the simulation ends when the 
Alfven waves have reached the upper boundary. 

Using the VAL C model we construct the temperature 
profile, and the density is determined by solving the hy- 
drostatic equilibrium equation. 

dP 

-QZ = ^^^^ 

The background equilibrium is shown in figure 4. 

All quanities were non-dimensionalised by using val- 
ues defined at z=0 from the VAL C model, thus 

r* = I50km (37) 

p* = 2.7 X 

V* = 65l5m/s 

f = Tis 



r = 6420/s: 

P* = 1.17xl0Va 
B* = lOOOG 

Conversion to variables with dimensions merely requires 
the multiplication of internal variables by these values. 

We model linear propagating Alfven waves by driving 
the horizontal velocity continuously at the bottom of 
the domain with a sinusoidal driving function, the ampU- 
tude of which did not affect the results. Typical driving 
velocities were 600 m/s and a range of 80 to 800 m/s was 
tested. 

As already mentioned in section 2 the Cowling resis- 
tivity is much larger than the Coulomb resistivity rj 
at chromospheric heights. In fact over the domain be- 
ing simulated, the calculated value of tj using VALC val- 
ues for the plasma never exceeds the value of numerical 
roundoff in the code, and hence all simulations are run 
with T] = 0. 

The dissipation of these waves can then be investigated 
by observing the change in amplitudes of the perturba- 
tions. Figure 5 shows a typical profile of the perpendic- 
ular velocity and magnetic field perturbations with height 
for a driving frequency of 0.07 Hz, for the case when no 
damping mechanism is applied and when coUisonal fric- 
tion damping is applied with the given profile of rjc (figure 
2). As in figure 4 the value of a is 0.3. 

By comparing the propagation of Alfven waves in the 
presence of a partially ionised region and without, we can 
study the efficiency of damping due to coUsional friction 
of ions and neutrals in the partially ionised chromosphere. 

As can be seen from figure 5, the gravitational strat- 
ification causes an increase in amplitude of the velocity 
pertubation and a decay in the corresponding magnetic 
field pertubations. This is a consequence of energy con- 
servation in the travelling wave. For a rigorous mathe- 
matical treatment of this see Moortel and Hood (2004). 
This makes a direct measuremnt of the efficiency of the 
coUisional frictional damping non-trivial. 

In order to obtain an estimate for the efficiency of col- 
lisional friction damping we compare the Poynting flux 
carried by the waves in the case when no damping mech- 
anism is present and when there is damping due to ion- 
neutral collisions. The time-averaged Poynting flux for 



7 



Figure 5: Velocity (top slides) and magnetic field (bot- 
tom slides) perturbations against height for a typical run 
with / = Q.QlHz, The left hand sUdes are for the case of 
r]c = and the right hand slides are for the rjc profile given 
by figure 2. The velocity and magnetic field are given in 
dimensionless units 



the Alfven waves is given by 



<5 >= 



^0 < b^v^ > 



(38) 



where Bq is the vertical background magnetic field and 
bx, Vx are the perturbations in magnetic field and velocity. 
By calculating the ratio of Poynting flux at a height above 
the partially ionised region of the model chromosphere 
we can estimate the efficiency of the damping mecha- 
nism. The ratio of damped to undamped Poynting flux 
at a height of 2500km above the surface (i.e. above the 
region of high Cowling conductivity) is given by 



< .V > 



dumped 



^ ^undamped 



(39) 



and for direct comparison with analyic results the effi- 
ciency of damping is estimated simply by 



= i-r. 



(40) 



Figures 6, 7, 8 and 9 show this estimate as a function of 
frequency, along with the estimate from the analytic ap- 
proach (solid lines). The four different plots are for the 
four different magnetic field profiles, with a being 0.3, 
0.35, 0.4 and 0.6 respectively. As the power law for the 
magnetic field changes so does our estimate for the damp- 
ing efficiency, as the value of rjc is dependant on |B|. 
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Figure 6: Estimates of the efficiency of damping due to 
ion-neutral collisions for the magnetic field profile with 
power law a = 0.3. The solid line is the estimate due 
to analytic approches and the triangles represent the esti- 
mates obtained from numerical data. 
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Figure 7: Estimates of the efficiency of damping due to 
ion-neutral collisions for the magnetic field profile with 
power law a = 0.35. The solid line is the estimate due 
to analytic approches and the triangles represent the esti- 
mates obtained from numerical data. 
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Figure 8: Estimates of the efficiency of damping due to 
ion-neutral collisions for the magnetic field profile with 
power law a = 0.4. The solid line is the estimate due 
to analytic approches and the triangles represent the esti- 
mates obtained from numerical data. 
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Figure 9: Estimates of the efficiency of damping due to 
ion-neutral collisions for the magnetic field profile with 
power law a = 0.6. The solid line is the estimate due 
to analytic approches and the triangles represent the esti- 
mates obtained from numerical data. 



The numerical data agrees best with the analytical es- 
timates when a is smallest (see figure 6). For this profile 
we can see that the two approaches agree in the middle of 
the frequency range, but disagree at low and high frequen- 
cies. The differences at low frequencies are due to limi- 
tations in the procedure outined in equations (38)-(40). 
The differences at high frequencies are due to the fact that 
the analytic approach is only applicable when the driv- 
ing frequency is well below the critical value (equation 
257), i.e the linear damping approximation. The differ- 
ence between numerical and analytic data increases as we 
increase the value of a (figures 6, 7, 8 and 9). By in- 
creasing this power in the magnetic field profile (equation 
25), we are decreasing the effect of ion-neutral collsional 
friction damping, as rjc is now smaller. This increases the 
relative importance of any other effects on the wave am- 
plitude. Thus the difference in numerical and analytic data 
is due to the fact that the numerical approach includes the 
stratification effects on the wave amplitudes, whereas the 
analytic approach does not. 

From both analytical and numerical data we can see 
that Alfven waves with frequencies below 0.01 Hz are 
unaffected by this damping mechanism, and propagate 
through the partially ionised region with little diffusion. 
Also waves with frequency above 0.6 Hz are completely 
damped by this mechanism. 

5 Conclusions 

We have estimated the efficiency of the partially ionised 
layer of the solar chromosphere in the damping of Alfven 
waves generated at the surface. The estimates were based 
on analytic and numerical approaches, which agreed in 
the linear damping approximation. The damping mecha- 
nism was the collisional friction between neutral and ion 
species in a model Hydrogen plasma with temperature 
and density values taken from the VAL C model of the 
quiet Sun ( IVernazza et al 19811 ). 

Alfven waves of frequencies above 0.6 Hz were com- 
pletely damped by the partially ionised layer in the chro- 
mosphere, whereas waves of frequency below 0.01 Hz 
were unaffected by the presence of neutrals and experi- 
enced no damping due to this mechanism. This lower 
result agrees with work conducted on MBP's and power 
spectra of horizontal motions at different heights in the 
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atmosphere (fCranmer and Ballegooijen in press]). They 
showed that using the WKB approximation the power 
spectra for motions with frequencies below 0.01 Hz are 
essentially unchanged as one progresses up the chromo- 
sphere to the transition region. This suggests that there 
is very little damping of low frequency waves due to any 
kind of dissipation mechanism in the solar chromosphere. 

Although higher frequency waves are difficult to ob- 
serve directly, photospheric motions can theoretically 
generate a large spectrum of Alfven waves. The fact that 
the partially ionised layer completely damps any waves 
above 0.6 Hz for our magnetic field models means that 
high frequency waves in the upper atmosphere must have 
been created by other sources than photospheric motions. 

This work is based on upward travelling waves gen- 
erated at the photospheric level. Downward travelling 
waves from the corona would be reflected at the den- 
sity contrast above the chromospheric region where the 
Cowling conductivity is large and it is unclear whether 
this damping mechanism would be important. The case 
of downward propagating waves will be subject to further 
investigation. 

The numerical simulations were performed in the non- 
ideal MHD approximation with an additional term relat- 
ing to the Cowling conductivity in the generalised Ohm's 
law (equation 16). Previous work on Alfven wave prop- 
agation in the lower solar atmosphere has used the WKB 
assumption, as used by (I DePontieu 19991 ). This assumes 
that the change in wavelength, as well as variables, is 
small over a typical wavelength. They estimated that in 
this regime the ion-neutral collisions in the chromosphere 
damped waves of frequencies above 0. 1 Hz. Our results 
are in broad agrement with the results of De Pontieu et 
al (2001) but as is clear from figures 7,8 and 9 the result 
is actually sensitive to the spreading out of the flux tube. 
Over the range of values of a, chosen to match estimates 
of |B| from observations, the damping is effective over a 
range 0. 1 to 0.6 depending on a. 

The model atmosphere was assumed to be magnetised, 
and thus the Generalised Ohm's law did not include the 
Hall term. In the upper photosphere, the electrons are 
tightly bound to the magnetic field whereas the ions are 
not. In this case there is a seperation electric field due to 
the neutrals drag on the non-magnetised ion, which should 



be taken into account in the one-fluid MHD equations. 
However, we showed that the damping of Alfven waves 
is efficient at heights of 1000km to 2000km and in this 
region the plasma can be regarded as magnetised. 

The form of the dissipation efficiency calculated from 
numerical data diff'ered from that obtained from linear an- 
alytic appoximations. Small differences could be seen at 
low frequencies (due to errors in the estimation of small 
damping decrements). Differences also occured at higher 
frequencies, although both analytic and numerical esti- 
mates of the efficiency must converge to 1 at high fre- 
quencies, the difference at higher frequencies is due to 
the fact that the linear damping approximation is not valid 
(equation 27). 
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